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1.  INTRODUCTION 


In  defense  preparedness  studies,  decisions  requiring  evaluation  of  the  effectiveness  of 
existing  and  conceptual  weapon  systems  often  arise.  Such  evaluations  usually  require 
some  measure  of  the  vulnerability  of  the  subject  systems  relative  to  a  specified  scenario 
and  mission  profile.  For  example,  if  one  is  considering  which  weapon  would  be  most 
effective  against  an  enemy  tank  of  the  1990's,  a  great  deal  of  uncertainty  surrounds  the 
description  of  this  future  threat.  Current  vulnerability  models  require  as  input  a 
detailed  computer  description  of  the  system  along  with  a  configuration  of  its  critical 
components. 

Present  methods  of  vulnerability  analysis  provide  a  point  estimate;  eg.,  a 
probability  of  kill  or  a  vulnerable  area,  with  no  corresponding  measure  of  uncertainty 
that  is  inherent  in  the  estimate.  Attempts  to  model  this  uncertainty  stochastically  have 
been  made  by  Dotterweich  1  and  Taylor  and  Bodt.  2 

In  th  is  document,  we  will  illustrate  how  vagueness  or  uncer’ainty  can  be  modi  led 
by  fuzzy  sets  and  how  this  uncertainty  is  reflected  in  the  final  measure  of  vulnerability. 
Our  emphasis  will  be  on  the  computation  of  vulnerable  area  of  a  surrogate  system  with 
uncertainty  present  in  the  cell  probabilities.  The  uncertainty  in  the  cell  probabilities 
can  be  attributed  to  uncertainty  in  assigning  component  kill  probabilities,  to  error  in 
engineering  judgement  in  determining  critical  areas  of  components  and  component  kill 
criteria  and  to  unknown  synergistic  effects.  Other  contributions  to  uncertainty  are  of 
course  present,  starting  from  the  target  description  through  a  myriad  of  subjective 
decisions  and  engineering  judgements  inherent  in  the  vulnerability  analysis. 

For  purpose  of  completeness  and  comparison,  a  point  estimate  of  the  vulnerable 
area  of  the  surrogate  system  will  also  be  determined  as  it  is  currently  done  in 
vulnerability  analysis.  A  point  and  interval  estimate  of  the  vulnerable  area  will  be 
determined  by  a  statistical  procedure  known  as  "bootstrapping."  Finally,  a  fuzzy 
vulnerable  area  will  be  calculated  and  the  different  approaches  compared. 


2.  DESCRIPTION  OF  A  VULNERABLE  .YREA  COMPUTATION 

The  models  currently  used  for  computing  vulnerable  area  of  armored  vehicle- 
require  a  number  of  detailed  inputs,  including  a  geometric  description  of  the  target 
specifying  the  location,  thickness,  and  obliquity  of  armor  plates  and  the  location  and 
size  of  major  components.  For  a  specified  direction  of  attack,  impact  locations  on  the 
target  are  selected  by  superimposing  a  rectangular  grid  and  choosing  one  point  of 
impact  within  each  grid  cell.  Then,  at  each  point  of  impact,  a  ray  is  traced  through  the 
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target  model  in  the  attack  direction  and  a  list  of  the  components,  barriers,  empty 
spaces,  and  other  structures  encountered  along  the  ray  is  compiled. 

For  each  component  encountered  along  the  ray,  a  fragment  mass  and  velocity 
together  with  an  estimate  of  fragment  shape  is  used  to  predict  a  hole  area  in  the 
component.  A  kill  criteria  in  the  form  of  a  hole  of  minimum  cross-sectional  area  in  a 
"sensitive  region"  of  the  component  is  required  in  order  to  regard  the  component  as 
inoperable.  If  the  predicted  hole  area  exceeds  the  minimum  hole  area,  then  a 
probability  of  kill  Pk  |h  j  is  assigned  to  the  jth  component  encountered  by  the  ith  ray.  In 
determining  probability  of  kill  given  a  hit  on  a  component,  the  vulnerability  analyst 
estimates  the  "presented  area"  of  the  component  and  the  presented  area  of  the  sensitive 
region  of  the  component.  The  estimated  probability  of  kill  is  taken  as  the  ratio  of  the 
presented  area  of  the  sensitive  region  to  the  presented  area  of  the  component  (An 
overall  component  kill  probability  is  obtained  by  summing  the  sensitive  area  from 
several  aspects  and  dividing  by  the  corresponding  total  component  presented  area.) 

If  several  components  are  encountered  along  a  ray,  and  determined  to  be  killed, 
then  the  probabilities  Pk|h  are  combined  under  an  independence  assumption  to 
calculate  a  probability  of  kill.  Pk|h,-  associated  with  the  ith  cell.  This  probability  is 
then  multiplied  by  the  area  of  the  cell.  A  j,  to  yield  a  "cell  vulnerable  area";  cell 
vulnerable  areas  are  then  summed  over  all  the  cells  in  the  grid  to  obtain  an  estimated 
target  vulnerable  area 

Av  =  EPkih.Ai  (2.1) 

i 

for  the  given  attack  aspect.  Hafer  and  Hafer  3  and  Nail  4  provide  further  details 
concerning  computation  of  vulnerable  area. 

3.  ESTIMATION  OF  WLNERABLE  .AREA  AND  PROBABILITY  OF  KILL 

Consider  the  conceptual  item  of  military  hardware  illustrated  in  Figure  1.  Without 
loss  of  generality,  the  component  edges  are  assumed  to  be  parallel  with  the 
superimposed  4x4  rectangular  grid.  The  numbers  in  the  ith  cell  represent  the 
probability  Pj,  j  h  that  the  target  will  be  killed  should  a  prescribed  round  of  ammunition 
fired  from  a  fixed  range  and  orientation  impact  in  that  location. 

The  vulnerable  area  Av  for  the  target  is  calculated  by 

\  S  r* It  I  h ,  Ai  =  A  v  Pk  I  h ,  (3.1) 

i  i 

^  T  F  Hafir  and  .4  S  Hafir,  ~  \  utnerabihly  Analyst  Jot  Surface  Target?  (VAST) ARBRL-TR-0S1S). 
Bathetic  R'narch  Laboratory  1979 

^  C  L  Sail,  *  Vulnerability  Analyst*  for  Surfart  Target*  (VAST),  Revition  Computer  Science * 
Corporation  Report  CSC-  TR-b £■  57^0.  !9bC 


where  i  is  indexed  over  (he  cells  of  the  grid.  In  t h is  example,  the  cell  areas  A,  are 
identical  and  equal  to  A.  The  overall  probability  of  kill  Pk  is 


where  Ap  is  the  total  presented  area  of  the  target.  For  the  conceptual  target  shown  in 
Figure  1,  the  exact  values  of  vulnerable  area  and  probability  of  kill  are  Av  =  105.1  and 
Pk  =  0  41 


In  practice,  the  values  Pk|h>  are  unknown,  and  estimates  Pk|h, 
obtain  the  approximations 


are  required  to 
(.13) 

(.3,1) 


Determining  valid  Pk  |  h  estimates  is  a  persistent  problem  because  of  the  difficulty  in 
modeling  the  underlying  damage  mechanisms,  and  is  an  additional  source  of 
uncertainty.  The  bootstrap  procedure  to  be  introduced  in  Section  4  assumes  the 
individual  Pk  |  h ’s  are  accurately  represented,  and  does  not  attempt  to  assign  a 
component  of  error  due  to  this  uncertainty.  This  problem  is  however  addressed  in 
See tit'ii  7  where  cell  Pkjhi’s  are  considered  to  be  fuzzy  numbers. 


4  THE  BOOTSTRAP 

The  bootstrap  is  a  technique  for  data  analysis  whose  goal  is  to  extract  information 
from  a  set  of  data  through  repeated  inspection.  An  expository  article  by  Diaconis  and 
Efron.  5  an  intermediate-level  paper  by  Efron  and  Gong.  6  or  a  more  advanced  paper  by 
Efron  '  serve  as  an  introduction  to  the  bootstrap. 

Suppose  we  want  to  estimate  some  attribute  9  of  a  population  X  =  { x , }  ,  £  j. 
where  the  index  set  I  may  be  finite  or  infinite.  To  obtain  an  estimate  9  of  6  ,  a  random 

n 

sample  {x.}  is  taken  from  X  and  a  statistic  0(x,.  •  -  ,  xn)  is  evaluated.  For  example. 

i-=  l 

if  9  =  /i  (the  population  mean),  then  estimates 

P  Diaconit  and  B  Efron,  Computer- Jntencive  Method/  «n  Stati/tic/.  Scirntific  American,  (1 !  C 
ISO 

^  B  Efron  and  G.  Gong,  A  Leisurely  Look  at  the  Bootstrap,  th e  Jackknife.  and  Cro//-\'ahdaticn.  Th ■ 
American  Statistician,  S7(19SS/,S6-4‘ 

^  B  Efron.  Bootstrap  Method/:  Another  Look  at  the  Jackknift,  Annalf  of  Statistic/.  7 i  1979/.1-C6 


arc  commonly  used;  if  6  =  <r2  (the  population  variance),  then 

n 

is  an  appropriate  statistic.  If  random  samples  {xj  are  repeatedly  drawn  from  X 

i«=  i 

and  a  statistic  0{  •  )  evaluated,  a  distinct  value  of  6(  •  )  may  he  determined  for  each 
sample  drawn.  These  values  are  from  the  so-called  sampling  distribution  of  6(  •  )  :  0(  ■  ) 
is  itself  a  random  variable  under  this  sampling  plan. 

Knowledge  of  the  sampling  distribution  of  0(  )  permits  questions  of  accuracy  to  be 

answered,  confidence  intervals  to  be  constructed  and  hypotheses  to  be  tested.  I'nder 
appropriate  assumptions,  and  for  some  functional  forms  of  6{  ■  )  ,  the  sampling 
distribution  may  be  determined  analytically;  however,  in  most  instances  it  cannot. 
Furthermore,  the  option  of  repeatedly  drawing  random  samples  to  construct  an 
empirical  approximation  to  the  sampling  distribution  for  6  is  usually  not  available. 

n 

Typically,  only  a  single  sample  of  data  { x , }  exists  for  analysis. 

The  essence  of  the  bootstrap  procedure  for  the  application  here  is  as  follows:  If 

n 

samples  with  replacement  from  {x, J  are  taken  to  generate  '"bootstrap"  samples 

i—  1 

{x,*}  and  6  =  0(x,* . xn*)  is  determined  a  large  number  of  times,  the 

1=1 

distribution  of  bootstrapped  6  provides  an  approximation  to  the  unknown  sampling 
distribution  of  6  .  (This  procedure  provides  a  Monte  Carlo  approximation  to  the 
distribution  of  the  nonparametric  maximum  likelihood  estimate  of  6  .) 

5.  BOOTSTRAP  DATA  ANALYSIS 

A  random  point  of  impact  and  its  corresponding  i  h  were  determined  for  each 

16 

cell  on  the  target  in  Figure  1.  These  sixteen  values  X=  {Pk|h}  form  the  data 

base  for  all  sub.-equent  analyses-  conventional,  bootstrap,  and  fuzzy.  In  this  instance, 
the  sixteen  values  formed  the  population  that  was  sampled  with  replacement  to 
determine  bootstrap  estimates  A*  of  the  vulnerable  area  A^.  Estimate  Av  was  based 
on  a  sample  equal  to  the  number  of  cells  in  the  overlaid  grid.  One  thousand  samples  of 
sixteen  were  drawn  from  X,  and  one  thousand  estimates  Ay  computed.  The 
approximation  to  the  distribution  of  Av  is  shown  in  Figure  2. 
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Fir.rr-.L  :  bodtsthaf  estimate  of  the  sampling  distriovtion  of  .4. 


The  histogram  in  Figure  2  is  a  bootstrap  estimate  of  the  sampling  distribution  of 
.  The  median  of  this  distribution,  which  we  will  take  as  an  estimate  of  is 
=  132.8  (compared  to  the  true  vulnerable  area  \  =  105.1).  An  approximate 
central  90T  confidence  interval  for  A^,  based  on  the  bootstrap  is  [99.3,  163.6]. 

In  this  example,  the  median  of  the  bootstrap  distribution  is  considerably  larger 
than  the  true  vulnerable  area.  This  is  because  the  sample  X  which  forms  the  basis  for  all 
the  bootstrap  calculations  estimates  A*.  to  be  133.2;  i.e.,  the  procedure  operates  on  a 
sample  providing  a  biased  estimate  of  Ay  .  The  bootstrap  has  no  way  to  compensate  for 
an  unknown  initial  bias,  and  in  this  instance  led  to  a  bootstrap  estimate  A^  =  132.8 
This  bias  has  no  effect  on  the  estimate  of  variance  of  A^.  since  the  error  is  in  location 
The  location  error  can  be  controlled  in  this  application  by  using  a  finer  grid  on  the 


6.  VAGUENESS  IN  YIXNERABIL1TY  ANALYSIS 


As  indicated  in  Section  2,  to  perforin  a  vulnerability  assessment  of  a  given  system, 
it  is  often  necessary  to  develop  an  elaborate  computer  description.  Information  about 
the  system  may  be  varied  in  form  and  in  abundance,  ranging  from  intelligence  data  to 
fully  documented  commercial  products.  The  system  may  be  foreign  or,  in  some 
instances,  conceptual.  One  of  the  most  subjective  parts  of  the  assessment  involves  what 
the  vulnerability  analyst  may  refer  to  as  "engineering  judgement."  Knowledge  of  a 
component’s  function,  for  example,  may  cause  the  analyst  to  predict  its  size  and 
location  in  a  system  for  which  little  information  is  available.  The  analyst  may  be  called 
upon  to  judge  armor  thickness,  material  and  thickness  of  component  structures,  and  to 
specify  kill  criteria  for  individual  components,  as  well  as  their  sensitive  regions.  The 
analyst  assesses  the  effect  of  different  types  and  degrees  of  damage  to  a  system  and  its 
components  and  infers  the  systems  "loss  of  function."  Given  damage  to  a  component  or 
subsystem,  the  analyst  may  infer  a  resulting  loss  of  mobility,  firepower,  etc.  The 
meaning  of  some  of  these  terms  and  appropriate  measures  for  their  quantification  are 
vague  or  imprecise;  even  vulnerability  analysts  differ  as  to  the  precise  meanings  and 
interpretations. 

At  many  places  in  a  vulnerability  analysis,  the  analyst  is  required  to  act  upon 
information  which  is  vague  or  imprecise  and,  in  some  cases,  where  information  is  not 
available  at  all.  Some  of  the  uncertainty  present  in  the  process  is  not  due  to 
randomness  and  therefore  probability  theory  may  not  provide  a  completely  satisfactory 
modeling  tool.  The  theory  of  fuzzy  sets  provides  a  rationale  for  dealing  with  sources  of 
uncertainty  which  are  non-random.  A  brief  introduction  to  fuzzy  sets  is  given  in 
Appendix  A.  In  the  following  section  we  illustrate  the  application  of  fuzzy  sets  to 
vulnerable  area  calculation. 


7.  CELL  PROBABILITIES  AS  FI  ZZY  NUMBERS 

The  vulnerable  area  model  MAST  requires  as  part  of  its  input  data  a  collection  of 
rays  corresponding  to  fragment  trajectories  for  given  angles  of  attack.  These  rays  are 
represented  by  lists  of  components,  armor  and  barriers  th3t  would  be  encountered  by 
the  fragments.  Along  each  ray  or  trajectory  the  mass  and  velocity  of  the  fragment  as  it 
exits  obstacles  it  encounters  is  computed.  For  an  initial  fragment  mass,  velocity  and 
shape  at  the  target  surface,  the  residual  mass  and  velocity  are  calculated  from  the  so 
called  THOR  8  equations.  As  the  fragment  encounters  a  particular  component,  a  hole 
area  is  estimated  in  that  component. 


Project  THOR,' The  Resistance  of  t'oriour  Metallic  Mat  trial*  to  Prrforation  by  Steel  /Va*r?;«  is.'.c 
Empirical  Relationships  for  Fragment  Residual  \'elocit y  and  Residual  Wright Projrct  THOR  T-  rh 
Rtpori  So.  47-  Ballistic  Analysis  Laboratory,  Inst.  for  Cooperative  Research,  The  John."  Hcpk.\  • 
Inivernty.  April  /19P1/. 


As  indicated  previously,  there  is  a  great  deal  of  vagueness  and  subjectivity  involved 
in  the  target  description,  material  specification  and  thickness,  component  I*k|h 
evaluation  and  kill  criteria.  Methods  currently  in  use  for  vulnerability  assessment, 
including  the  bootstrap  procedure  discussed  in  Section  4,  do  not  take  these  uncertainties 
into  account.  Assumptions  regarding  the  fragment  shape,  barrier  material  and  thickness 
dictate  that  the  hole  estimate  is  actually  a  fuzzy  number. 

Incorporating  fuzziness  into  both  the  kill  criteria  and  the  component  conditional 
kill  probability  results  in  the  following  fuzzy  implication:  ’’’If  the  cross-sectional  hole 
area  in  the  component  is  about  £  cm2  or  greater,  then  the  component  conditional  kill 
probability  Pk  [  h is  B,"  where  B  is  now  a  fuzzy  number.  Thus,  we  have  a  fuzzy  hole 
area  estimate  and  a  fuzzy  implication  relating  kill  criteria  with  kill  probability.  If  these 
wore  all  ’’crisp,"  as  assumed  in  the  current  method  for  vulnerability  assessment,  one 
nrod  only  compare  the  predicted  hole  area  with  that  of  the  kill  criteria.  If  the  predicted 
hole  area  is  greater  than  the  required  value,  then  a  probability  of  kill  is  assigned.  This  is 
described  by  the  syllogism 

If  h  is  A,  then  p  is  B 
h  is  A _ 

p  is  B 

However,  when  A  and  B  are  fuzzy  and  we  have  h  is  A*,  then  we  must  use  the 
compositional  rule  of  inference  or  the  generalized  modus  ponens  of  fuzzy  logic  9  to  obtain 
the  conclusion.  This  can  be  represented  by  the  syllogism 

If  h  is  A,  then  p  is  B 
h  is  A* _ 

p  is  B* 

The  resulting  fuzzy  component  conditional  kill  probabilities,  when  combined,  result  in  a 
fuzzy  cell  probability  and  finally  a  fuzzy  vulnerable  area  estimate,  reflecting  the 
uncertainty  in  the  process. 

Zadch  10  defined  "If  A  then  B"  as  a  special  ease  of  the  fuzzy  implication  "If  A  then 
B  else  C.”  Notaf ionally,  we  represent  If  A  then  B  else  ('  as  (  AxB  )  (j  (  -AxC'  )  where 
indicates  negation  and  X  is  the  cartesian  product  of  the  two  fuzzy  sets.  If  (’  is  ihe 
empty  set,  then  If  A  then  B  is  written  (  AxB  ).  Recalling  that  the  cartesian  product  of 
two  fuzzy  sets  A  and  B  is  given  by 

/*a*b  (x-  >')  =  rnin{/<A  (x)-  f'B  (>')}•  ("  D 

9 

L  A.  Zadeh,  '  Outline  of  a  ,VfU-  Approach  to  the  Analyst?  of  Complex  System?  and  Deeirion  Process- s. 
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we  ran  derive  a  relation  between  the  hole  area  and  the  kill  probability,  namely 

Pr  (h,  P)  =  min{ji„  (h),  //P  (p)}  (7.2) 

•J 

whore  R  =  II  X  Pk(h  j. 

Example :  If  II  =  0. l/hj  +  0.5/h2  +  0.7/h3  +  1/h,  where  h;  <  hj  if  i  <  j,  and 
I\  |h „  =  0  I/P,  +  0.5/p2  +  l/l»3  +  0.5/p,  +  0.1/p5  ,  then  /iR  (h,  p)  is  given  by 


Pn  (h.  p) 

p 

Pi 

P2 

P3 

Pi 

Ps 

h 

h, 

0.1 

0.1 

0.1 

0.1 

0  1 

h2 

0  1 

0.5 

0.5 

0.5 

0.1 

^3 

0.1 

05 

0.7 

0.5 

0.1 

0.1 

0.5 

1.0 

0.5 

0.1 

Furthermore,  if  we  assume  that  the  hole  estimate  H,  is  given  by 

II,  =  0  3/h ,  +  l/h2  +  0.5/b3  +  0.1/h, 

then  we  would  like  to  combine  this  estimate  with  that  given  by  the  fuzzy  implication  to 
infer  something  about  Pk  |  h T°  proceed  further,  we  need  to  invoke  the  compositional 
rule  of  inference. 

Given  universes  l'  and  V,  a  fuzzy  set  A  on  l1  and  a  fuzzy  relation  R  on  I'xW 
characterized  by  //A  (u)  and  //R  (u,  v)  respectively,  the  compositional  rule  of  inference 
states  that  the  result  will  be  a  fuzzy  set  B  on  V  defined  by 

pb  (v)  =  u !  i‘.\  <u)  n  pr  (u- v) )  (“•'*; 

u 


and  written  as  B  =  A  o  R. 


.  * 


Continuing  the  above  example,  we  Lave 

H,  o  R  =  B 

with  D  denoting  the  fuzzy  set  resulting  from  this  composition.  Formally, 

Ph  (Pj)  =  U{  /'ll,  00  D  PR  (*V  Pi)}  (71) 

=  u  {0.3 n o.i,  i.ono.i.o.5no.i,o.i  no.i) 

=  y  {0.1.  0.1,  0.1,  0.1} 

=  0.1. 

Similarly.  //R  (p2)  =  (p3)  =  //q  ( p4)  =  0.5,  and  //R  (p5)  =  0.1;  hence  11,  o  R  =  B. 

If  11,  =  H,  and  both  are  non-fuzzy,  then 

HoR  =  Pk(h)  (7.5) 

in  agreement  with  the  classical  modus  ponens.  With  II,  ^  H  and  fuzzy,  v.e  h .  vc  a 
fuzzy  version  of  the  classical  modus  ponens  with  the  important  dilTerence  that  11,  and  11 
can  be  different.  Thus,  in  the  above  example,  we  have  a  hole  estimate  which  C 
different,  in  general,  from  the  required  kill  criterion. 

Similar  results  can  be  obtained  for  each  component  encountered  along  the  fragment 
trajectory.  These  kill  probabilities  must  be  ultimately  aggregated  into  a  cell  kill 
probability.  As  mentioned  in  Section  2,  current  methods  include  combining  component 
kills  as  independent  events;  however,  the  problem  of  aggregation  of  these  probabilit ins 
remains  open.  Some  methods  of  aggregating  fuzzy  sets  such  as  union  and  intersection 
are  presented  in  Appendix  A.  Zimmerman  and  Zysno  11  have  investigated  several  other 
methods.  For  vulnerable  area  calculations,  aggregation  of  the  fuzzy  sets  requires  further 
research;  therefore,  in  the  following  we  shall  assume  that  some  form  of  aggregation  has 
been  done,  and  the  resulting  probability  of  kill  associated  with  each  cell,  Pk  |  h  .  is  a 
fuzzv  number. 


5$ 
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1  H-J  Zimmerman  and  P.  Zysno,  " Latent  Connectives  in  Hitman  Decision  Staking,'  Fu::y  Sets  at :d 
System*  4(1956),  57-57. 
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8.  FI  ZZY  ARITHMETIC 

In  (his  section,  wo  consider  P^j^  *°  be  »  fuzzy  number  and  A  (o  lie  a  ori-p 
number;  to  simplify  (he  notation  we  will  denote  P;  =  Pk  |  (l  Equation  (3.1)  will  th<n  I .*■ 
evaluated  by  fuzzy  arithmetic  to  obtain  a  fuzzy  number 

AV  =  AYIV  (8.1) 

i 

Let  /'r,(x)  and  /ip.  (y)  be  normalized  membership  functions  of  I*,  and  P2  with 
support  in  [0.1],  The  definition  of  the  sum  of  two  fuzzy  numbers  P,  +  J\  is  pi w*n  by 


f'tr.  +  r)(,1=  max  (  min  (/tpi  (x),  /i  (y ))  j 
x  +  \  =  i  L  J 


I  V> 


-I 


\Y  e  will  choose  our  class  of  membership  functions  to  be  continuous  and  symmetric  about 
a  point  x0  such  that 


f‘P  (*)  = 


monotone  nondecreasing  x  <  x0 
monotone  nonincreasing  x  >  x0  . 


(83) 


In  particular,  consider  the  following  two  functions. 


I‘  p,(x)  = 


0 

f(x) 

1 


X  <  X, 

x,  <  x  <  x2 
x2  <  X  <  x3 


(8  1) 


and 


f‘pjy)  = 


0  y  <  y, 

g(y)  >  i  <  y  <  y2 
i  >'2  <  >  <  >3 


(8  3) 
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where  the  functions  ft  (x)  and  ft  (y)  are  symmetric  about  x3  and  y3,  respectively 
Dubois  and  Trade  show  how  to  operate  on  fuzzy  numbers  of  this  form  in  an  efficient 
manner.  In  Appendix  B  a  geometric  interpretation  useful  in  evaluating  (8.2)  for  our 
specific  choice  of  membership  function  is  detailed. 


The  product  of  two  fuzzy  numbers  A  and  P  is  defined  as 

p(t)=  max  [min  (//A  (x),  ft  (y ))  1 
xy  =  t  1  J 


(8  0) 
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In  this  instance,  A  is  considered  a  crisp  number,  which  has  as  its  membership  function 


''a(*>={  o  otherwise.  (X7> 

Since,  for  a  given  t.  (\)  =  0  for  x  ^  A,  the  only  value  under  consideration  is  x  =  A 
This  implies 

/'a  p(l)=  ™ax  [min  (1,  /ip  (>  ))]  (8.S) 

Ay  =  t  !■  ■* 

8.  FIZZY  DATA  ANALYSIS 

IS 

We  will  now  "fuzzify"  the  data  base  <Pk|  h  }  introduced  in  Section  5  for  the 

bootstrap  procedure  and  calculate  the  vulnerable  area  A^  using  fuzzy  arithmetic.  In 
choosing  the  membership  functions  /ip  ,  i=l,  •  •  •  ,16,  we  make  use  of  the  fact  that  the 
variance  of  the  estimate  Pkjh  is  Pk  j  h ( 1  -  Pk  j  h  )  although  no  statistical  properties  will 
be  invoked  in  the  ensuing  computation. 

Let  A,  =  1/2  Pk  |  h,  (1  -  Pk|  h,)  •  The  membership  function  for  P;  was  taken  to  be 

0  x  <  pk|h,  -  2Aj 

L,(x)  Pk,h.-2Ai<x<Pk|hi-Ai  (9.1) 

1  |  h,  _  <  x  <  I  h, 


where  L,(x)  is  the  line  determined  by  the  points  (Pk  |  h  -  2Aj,  0)  and  (Pk|h  -  A;.  1)  : 
also,  recall  that  //p ( x )  is  symmetric  about  Pk | h  . 


For  example,  the  fuzzification  of  PM  =  0.6  has  as  its  membership  function 


0  x  <  0.36 

25x/3  -  3  0.36  <  x  <  0.48  . 
1  0 .48  <  x  <  0.60 


(9.2) 
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From  Dubois  and  Prade  or  Appendix  B  the  fuzzy  Dumber  P  =  VPj  is  given  by 

i-l 


the  membership  function 


0  x<a 

/^p(x)  =  l(x)  a<x<b 
1  b<x<c 


(0.3) 


and  l(x)  is  a  line  segment  determined  by  the  points  (a,0)  and  (b,  1 ). 

Finally,  A*,  the  product  of  the  crisp  number  A  and  the  fuzzy  number  P  has  as  its 
membership  function 

0  x  <  Aa 

/iAJx)  =  L(x)  Aa  <  x  <  Ab  ,  (9.4) 

1  Ab  <  x  <  Ac 

where  L(x)  =  l(x/A)  on  (Aa,  Ab).  The  membership  function  for  A,  for  the  example  in 
Section  3  is  given  by 

0  X  <  93.7 

PaJx)  =  05lx  -  4.78  93.7  <  x  <  113.4  .  (9.5) 

1  113.4  <  x  <  133.2 

as  illustrated  in  Figure  3. 

Since  the  function  /iAjx)  is  symmetric  about  x  =  133.2,  values  of  vulnerable  area 
in  the  interval  (113.4,  153.0),  where  the  membership  function  takes  on  value  one,  are 
equally  acceptable.  The  vulnerable  area  obtained  from  conventional  calculation  on  the 
sixteen  random  shotlines  is  at  the  midpoint  of  this  interval,  while  the  true  vulnerable 
area  for  this  example,  \  =  105.1,  has  from  (9.5)  a  membership  value  of  0.58  . 

Since  the  fuzzy  vulnerable  area  depends  on  the  cell  probabilities,  the  fuzzy 
vulnerable  area  is  dependent  on  the  grid  size.  Grid  size  impacts  the  bootstrap  and 
conventional  calculations  as  well.  In  this  example,  if  the  number  of  cells  is  increased  by 
a  factor  of  four,  the  interval  of  values  for  which  /*A(x)  =  1  has  as  its  midpoint  the  true 
value. 


Furthermore,  the  interval  of  values  for  which  the  membership  function  equals  one 
depends  upon  the  interval  of  values  for  which  the  individual  celi  probabilities  have 
membership  values  equal  to  one.  We  assumed  in  (9.1)  a  specific  form  for  the  fuzzy  cell 
probabilities.  The  vulnerability  analyst  may  be  able  to  sharpen  these  bounds,  and  in  so 
doing,  sharpen  the  fuzzy  vulnerability  estimates. 


10.  SIWMARY 

In  this  report,  we  have  attempted  to  make  the  point  that  at  many  places  in 
vulnerability  analysis  there  is  a  margin  for  error  that  comes  about,  not  from 
randomness,  but  from  a  basic  inability  to  model  uncertainty.  We  have  considered,  for  a 
specific  example,  two  new  approaches  for  vulnerability  analysis:  (i)  a  fuzzy  set  approach 
and  (ii)  a  statistical  resampling  plan  known  as  the  bootstrap. 

The  membership  function  /^(x)  in  (9.5)  for  the  fuzzy  vulnerable  area  calculation  is 
illustrated  in  Figure  3.  The  fu;:y  procedure  seeks  to  quantify  the  uncertainty  inherent  in 
a  vulnerable  area  calculation  due  to  the  uncertainty  in  the  ce//Pk|h  estimates.  This 
uncertainty  stems  from  the  assignment  of  fuzzy  kill  probabilities  to  the  individual 
components  and  from  the  method  of  aggregation  into  a  cell  probability  of  kill.  It  has 
the  distinct  advantage,  however,  of  clearly  delineating  the  inherent  subjectivity. 

The  point  estimate  of  \  provided  by  the  current  method  of  vulnerability  analysis. 

=  133.2,  is  indicated  in  Figure  3.  The  fuzzy  vulnerable  area  is  symmetric  about  Av. 
but  considers  values  in  the  interval  (1 13.4,  153.0),  corresponding  to  values  of  A*  for 
which  the  membership  function  equals  one,  to  be  equally  acceptable.  This  interval  has 
the  appearance  of  a  statistical  confidence  interval,  but  does  not  convey  the  same 
information,  since  it  has  no  corresponding  probability  level  attached. 

As  discussed  in  detail  in  Section  5,  the  sampling  distribution  of  \  provided  by  the 
bootstrap  has  a  median  of  132.8  and  an  approximate  central  90^  confidence  interval  of 
[99  3.  1G3.6]  The  bootstrap  procedure  seeks  to  quantify  uncertainty  inherent  in  a 
vulnerable  area  calculation  due  to  the  variability  in  the  sampling  procedure  responsible 
for  selection  of  the  cell  Pk  j  h  estimates.  The  bootstrap  takes  the  cell  Pk  |  h  as  gi\ en.  and 
makes  no  attempt  to  model  uncertainty  in  the  values  themselves. 


FIGURE  S  MEMBERSHIP  FUNCTION  FOR  FUZZY  A 


The  estimate  is  subject  to  error  due  to  the  uncertainty  in  the  cell  I\  1 1 
estimate^.  the  variability  in  the  sampling  scheme  that  selects  the  cell  I\  |  }l(  estimate*, 
the  mesh  size  of  the  overlaid  target  grid,  and  likely  additional  sources.  Tim 
implementation  of  this  observation  would  suggest  bootstrapping  on  fuzzy  numbers.  The 
theoretical  development  is  not  sufficient  at  this  time  to  allow  this,  but  it  is  a  logical 
research  extension.  Both  the  fuzzy  set  approach  and  the  bootstrap  are  worthy  of 
implementation  into  the  vulnerability  estimation  regime.  Both  procedures  enhance  the 
vulnerability  estimates  as  they  are  now  provided  without  undermining  the  lore  of 
vulnerability  analysis  as  currently  conceived. 

In  Section  7  the  use  of  fuzzy  sets  in  expressing  the  uncertainty  in  the  hole  area 
required  in  a  sensitive  area  of  a  component  to  achieve  a  component  kill  was  illustrated 
The  fuzzy  implication  formed  by  the  fuzzy  criterion  and  the  fuzzy  component  kill 
probability  was  defined  by  the  cartesian  product  of  the  two  fuzzy  sets  (7.1).  Ollier 

definitions  exi*t  for  implication  and  compositional  rules;  eg.  the  recent  work  of 

1 1 

Mizumoto. 


/  *? 

Mi  :umoto,  'Fu::y  Rtaforung  Cndrr  ,Vrtr  Compcfitional  liultr  of  Infi  r< nn kybornrtr r,  l.'i 
107-117. 
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Zadeh.  14  noting  tint  ninny  classes  of  objects  encountered  in  t lie  real  world  nr.'  n  • 
precisely  defined,  proposed  the  concept  of  a  fuzzy  set  to  model  sijoJi  classes.  Tin  < )  ,  f 

all  real  numliers  which  are  "pinch  greater  than  one"  or  the  class  of  "severely  d.:in  eo  j 
vehicles,"  for  example,  do  not  constitute  classes  in  the  usual  mathematical  sense.  |  e//\ 
set  methodology  hold'  promise  for  vulnerability  analysis  and  damage  assessment  -  -m  • 
it  proiides  a  mechanism  by  which  subjective  information  ran  be  quantified  and  op.  rat.  I 
on.  with  results  reflecting  the  degree  of  imprecision  in  input  information  and  data  In  a 
closely  related  field,  Blockley  1  ,  lf’  and  Yao  1  ‘  have  successfully  applied  fuzzy  i<, 

problems  in  ci\il  engineering  and  structural  damage. 

The  fundamental  notion  of  a  fuzzy  set  i'  the  following: 

If  V  (the  universe  of  discourse)  is  a  collection  of  objects,  then  a  fu;:y  set  A  i-  a  “•!  of 
ordered  pairs 

A  =  {(X.  #*A  (X)  )  |  X  <  V  }  . 

where  // A  (\)  ,  called  a  membership  function,  quantifies  the  grade  of  membership  of  \  in 
A.  The  function  ft A  (x)  maps  x  into  a  membership  space  M  which  is  usually  taken  to  |." 
the  unit  interval  [0.1], 

Example.  1‘  =  Authors  of  this  report:  A  =  The  set  of  bald  men 

A  =  {(Malcolm.  0.1).  (Palmer.  0  7).  (Ralph.  0.0)} 

For  "crisp"  sets,  i.e.,  ordinary  sets,  the  membership  function  reduces  to  th" 
characteristic  function 

»x»  _  (  1  x  c  A 
/  A  |X|  -  }  0,  otherwise  ' 

Fuzzy  sets  need  not  be  finite  sets;  for  example,  in  vulnerability  studies  vehicle'  take  <>n 
"loss  of  function"  values  in  the  unit  interval  [O.lj.  One  possibility  for  a  fuzzy  >.t 
characterizing  "high  loss  of  function"  is  illustrated  in  Figure  1 


^ L  A  Zadth.  '  F uzzy  Site,’  Information  and  Control.  ill3Cf<i.  SSSSCS 

I  RtockU  y.  '  Tht  Half  of  Fuzzy  St  /.»  in  Civil  Engint  r  ring'  Fuzzy  St  t  and  Syrtcme.  Cl  19  79 
1  D  1  Dlc-klty  Tht  Satur •  of  Structural  l)<  sign  and  Softly.  Fill?  lloruoed  I  iriitr  d.  tl9<0i. 
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FIGl'RE  4  MEMBERSHIP  FUNCTION  FOR  'HIGH  LO£S-OF-Fl’NCTiON  • 


We  now  give  several  definitions  for  properties  of  fuzzy  sets  which  are  obvious 
extensions  of  corresponding  definitions  for  crisp  sets. 

A  fuzzy  set  A  is  empty  if  and  only  if 

fiA  (x)  —  0  for  all  x  t  U  . 

Two  fuzzy  sets  A  and  B  are  equal  if  and  only  if 
Ha  (x)  =  Hb  (x)  toT  *11  X  t  U  . 

The  complement  of  a  fuzzy  set  A,  denoted  by  Af,  is  defined  by 
HA‘  (x)  =  1  *  Ha  (x)  »H  x  «  U  . 

The  union  of  two  fuzzy  sets  A  and  B,  with  membership  functions  / ja  and 
is  a  fuzzy  set  C  whose  membership  function  is  given  by 

pc  (x)  =  max  {  pA  (x),  pB  (x)  },  xtU 

alternatively  written 

He  (x)  =  Ha  (x)  U  Hb  (x)>  x  i  U 

where  u  denotes  the  maximum  operator. 


Tho  intersection  of  two  fuzzy  sets  A  and  Bisa  fuzzy  set  C  given  by 
Pc  (x)  =  Pa  (x)  n  PD  (x),  X  i  U 
where  p)  denotes  the  minimum  operator. 

The  support  of  a  fuzzy  set  A  is  a  crisp  set  S(A),  given  by 
S(A)  =  {  x  |  /iA  (x)  >  0}  . 

A  fuzzy  set  A  is  said  to  be  normal  if  and  only  if 

sup  Pa  (x)  =  1 
x 

The  cartesian  product  AxB  of  fuzzy  sets  A  and  B  from  universes  U  and  V, 
respectively,  is  defined  by 

PajB  (x.  y)  =  min  {Pa  (x),  Pb  (y)}.  x  (  U,  y  <  V 

A  fuzzy  number  A  is  a  fuzzy  subset  of  the  real  line  R. 

A  fuzzy  number  A  is  said  to  be  convex,  if  for  any  real  numbers  x,  y,  z  t  H 
with  x  <  y  <  z 

Pa  (>’ )  >  Pa  WHPaW 

An  a-let  el  set  of  a  fuzzy  number  A  is  a  crisp  set  A0  defined  as 
A0  =  {x  I  fJA  (x)  >  a},  0  <  a  <  1  . 

The  fuzzy  set  A  is  commonly  written,  when  U  is  finite,  as 

A  =  PA  (x,)  /  x,  +  •  •  •  +  fiA  (xn)  /  xB  =  £  nA  (xj  /  Xj 

i 

where  V  denotes  the  union  of  fuzzy  sets;  and  when  U  is  not  finite  as 

A  —  /  pA  (x)  /  x  . 


The  algebraic  properties  of  fuzzy  numbers  under  the  four  arithmetic  operation^  of 
fuzzy  addition,  subtraction,  multiplication  and  division  have  been  extensively  stndi  d  by 
Mizumoto  and  Tanaka  18  and  Dubois  and  Trade.  1 J  For  example,  if  A  and  B  are  fuzzy 
numbers  with  membership  functions  pA  and  /iy,  then  the  sum  of  A  and  B  is  given  by 

/'A+n  (z)  =  U  O'aMP/'u  (>)) 
x  +  y  =  z 

=  U  </'a  (*)  fl  /' n  <z  -  x) ) 

x 

(Recall  that  (j  and  p|  represent  the  maximum  and  minimum  operators.) 

Example  :  If  l  =  {0,  1,  2,  ...  },  A  =  0.2/1  +  1/2  +  0.3/3  and 
B  =  0.1/2  +  1/3  4-  0  2/1.  then  z  can  take  the  values  3,  4,  5,  G  and  7  and 

/'a+d  H)  =  max  {//A  (1)  n  /'u  (3)-  /‘a  (2)  Pi  f'D  (2>} 

=  max  {0.2  n  I,  1  p  0.1}  =  ntax  {0.2,  0.1} 

=  0.2. 


Fuzzy  addition  is  clearly  more  cumbersome  than  crisp  addition;  liovxur. 
summation  may  be  readily  implemented  on  the  computer. 


1 1 

'M  Mizumoto  and  l\  Tanaka.  " Some  Properhe s  of  Fuzzy  .\’u»»i*r  r.«."  Advance.-  in  Fuzzy  S-t  77. ■  . 
Application ?,  A/.  A/.  Gupta.  /!.  h.  Ragadt,  R.  R  Yager.  Editor?.  Sorth-HAir.nd  (I97J). 

P D .  Dubois  and  II  Pradi,  op.  fit. 

CD  ■ 

D.  Dubois  and  //  f’radt,  Fuzzy  Fit?  and  Systems:  Theory  and  Application*.  Arcd>r:ic  Press.  IIP', 


Appendix  B.  A  GEOMETRIC  INTERPRETATION  OF  FUZZY  ADDITION 


In  order  to  determine  the  sum  of  Pj  and  P2,  we  will  use  a  geometric  interpretation 
to  evaluate  (8.2).  On  an  xy-coordinate  system,  graph  pPj  (x);  for  the  membership 
function  /ipj  (y),  graph  the  independent  variable  along  the  y-axis  as  shown  in  Figure  5. 
For  a  given  t,  consider  the  line  x  +  y  =  t;  the  values  that  satisfy  the  equation  of  this 
line  will  be  used  to  evaluate  /i^pj  +  pjj  (t). 


FIGURE  b  MEMBERSHIP  FUNCTIONS  FOR  FUZZY  P,  AND  Pa 

If  t  <  x,  +  y,  (Figure  6(a)),  then  either  ppi  (x)  =  0  or  ppj  (y)  =  0;  this  implies 

min  [pPl(x)  ,  Ppa(y)j  =  0.  (B.I) 


From  (8.2),  p(pi  +  Psl  (t)  =  0  .  By  symmetry,  an  analogous  result  holds  for  t  >  xb  +  y5. 

Now  consider  t  in  the  interval  x2  +  y2  <  t  <  x4  +  y4.  Then  the  line  x  +  y  =  t 
passes  through  the  rectangle  determined  by  the  line  segments  x2  x4  and  y2  v4  (Figure 
6(b)).  At  every  point  (x,y)  in  the  rectangle,  and  in  particular,  at  every  point  on  the  line 
x  +  y  =  t  within  the  rectangle,  the  value  of  both  membership  functions  is  one; 
therefore, 

P(Pl  +  ps)  W  —  1  for  x2  +  y2  <  t  <  x4  +  y4  .  (B.2) 

Finally,  consider  t  in  the  interval  Xj  +  y,  <  t  <  x2  +  y2.  On  the  line  x  +  y  =  t 
(Figure  6(c)),  as  x  takes  on  values  zero  to  x2,  /Jpj(x)  is  nondecreasing  from  zero  to  one; 
simultaneously,  the  values  of  #ipj(y),  y  =  t  -  x,  are  nonincreasing.  Thus, 

29 


■  1,  » - *- 


T 


until  //pi(x)  =  /ipj( y),  then 


min  ^Pl(x),  /<Ps(y)J  =  #ipj(y) 


Therefore, 


/i(r.4rj(t)  =  #ipJx‘(t)). 
where  x*  (t)  satisfies  the  equation 

<*•>-#«*  <t  -  o . 

A  similar  argument  holds  for  x4  +  y4  <  t  <  x5  +  ys. 

Suppose  f  (x)  and  g  (y)  have  the  form  f  (x)  =  ajX  +  bj  and  g  (y) 
Then  the  x*  that  satisfies  f  (x*)  =  g  (t  -  x*)  is  given  by 

.  *2  .  .  ^2  "  bj 


=  a2y  +  b2. 


x*  = 


- 1  +  - 

aj  +  ao  i]  *f  &2 


Thus,  the  sum  of  two  fuzzy  numbers  from  the  class  (8.4-5)  maintains  the  same  analytic 
structure;  that  is,  line  segments  are  invariant.  However,  the  support  is  now  contained  in 
[0,2].  This  argument  extends  naturally  to  the  sum  of  n  fuzzy  numbers. 


<  *,*7, 


FIGURE  «(»)  CASE  I  I  <  +  V, 


■t 


V  >  *j4T, 

( <  *4**4 


I  *1**  *4*1 

FIGURE  6(b)  CASE  2  %  +  Vs  <*<**+  V* 


FIGURE  6(c)  CASE  3 
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